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Abstract 

Background: Understanding how biodiversity is shaped through time is a fundamental question in biology. Even 
though tropical rain forests (TRF) represent the most diverse terrestrial biomes on the planet, the timing, location 
and mechanisms of their diversification remain poorly understood. Molecular phylogenies are valuable tools for 
exploring these issues, but to date most studies have focused only on recent time scales, which minimises their 
explanatory potential. In order to provide a long-term view of TRF diversification, we constructed the first complete 
genus-level dated phylogeny of a largely TRF-restricted plant family with a known history dating back to the 
Cretaceous. Palms (Arecaceae/Palmae) are one of the most characteristic and ecologically important components 
of TRF worldwide, and represent a model group for the investigation of TRF evolution. 

Results: We provide evidence that diversification of extant lineages of palms started during the mid-Cretaceous 
period about 100 million years ago. Ancestral biome and area reconstructions for the whole family strongly 
support the hypothesis that palms diversified in a TRF-like environment at northern latitudes. Finally, our results 
suggest that palms conform to a constant diversification model (the 'museum' model or Yule process), at least 
until the Neogene, with no evidence for any change in diversification rates even through the Cretaceous/ 
Paleogene mass extinction event. 

Conclusions: Because palms are restricted to TRF and assuming biome conservatism over time, our results suggest 
the presence of a TRF-like biome in the mid-Cretaceous period of Laurasia, consistent with controversial fossil 
evidence of the earliest TRF. Throughout its history, the TRF biome is thought to have been highly dynamic and to 
have fluctuated greatly in extent, but it has persisted even during climatically unfavourable periods. This may have 
allowed old lineages to survive and contribute to the steady accumulation of diversity over time. In contrast to 
other plant studies, our results suggest that ancient and steady evolutionary processes dating back to the mid- 
Cretaceous period can contribute, at least in part, to present day species richness in TRF. 



Background 

Tropical rain forests (TRF) are the most biodiverse ter- 
restrial ecosystems on the planet [1]. They are charac- 
terised by a closed, multilayered canopy dominated by 
flowering plants (angiosperms [1]) and occur only in 
frost-free areas with high mean monthly temperatures 
and precipitation, and low seasonality [2]. Today, TRF 
covers just 7% of the Earth's surface [3] in equatorial 
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zones of the Americas, Africa and the Indo-Pacific, and 
is highly threatened by human activity [4]. The origin 
and evolution of species-rich biomes raises fundamental 
questions in evolutionary biology [5] and, as such, the 
diversification of TRF has been much debated [3,6] . 

Even though it is generally agreed that TRF is a relatively 
old biome, the location and timing of its origin remain 
uncertain mainly because the fossil record for tropical 
regions is highly incomplete, especially during the Cretac- 
eous [7,8]. Some direct [9] and indirect [10] evidence sug- 
gests that TRF was present in the mid-Cretaceous period 
(100 million years ago (Ma)) at middle paleolatitudes (for 
example, Laurasia) while other studies indicate that the 
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first paleofiora attributable to modern day TRF are found 
in the Early Paleocene of North America [11] and Late 
Paleocene of South America [7,8,12] and Africa [13-15]. 

Whereas previous views suggested that the TRF biome 
has been ecologically stable over long periods of time 
[16], more recent data indicate that it is highly dynamic 
[17] having fluctuated both in extent [14,18] and in the 
diversity of plants that it sustains [7,8]. These views 
have led to three general evolutionary hypotheses that 
explain the high levels of present day species diversity 
found within TRF: (i) early, rapid speciation in response 
to favourable climatic conditions followed by a decelera- 
tion of diversification rates due to global cooling [19] 
and contraction of TRF (referred to here as the 'ancient 
cradle model', see [10,20,21]); (ii) constant diversification 
rates coupled with low extinction rates leading to a gra- 
dual accumulation of lineages in response to a long-last- 
ing and stable tropical ecosystem (the 'museum model', 
see [16,22]); and (Hi) an increase in diversification rates 
towards the present in response to climatic, tectonic or 
biotic changes (the 'recent cradle model', see [23-25]). 
These hypothetical processes result in alternative lineage 
accumulation through time and thus different patterns 
of inferred branch length distributions (Figure 1). Sev- 
eral phylogenetic studies of TRF plant groups have pro- 
vided evidence in support of the recent cradle model of 
diversification [25-27]. However, these studies were 
restricted to low taxonomic levels (species) and thus do 
not enhance our understanding of how these hypotheses 




Time since origin 

Figure 1 Lineage-through-time (LTT) plots for three alternative 
hypothetical diversification models of tropical rain forest (TRF) 
evolution. (1) Decrease in diversification rates since origin with 
early radiation; 'ancient cradle model'; (2) constant diversification 
rate, 'museum model'; (3) increase in diversification rates since origin 
with recent radiation; 'recent cradle model'. 



might apply throughout the entire history of TRF, for 
example, on the long-term diversification dynamics of 
TRF. In contrast, studies of the early diversification of 
TRF plant lineages that permit tests of the above 
hypotheses are rare. 

In the absence of a complete fossil record for tropical 
rain forests, family-level diversification analyses of large 
pantropical angiosperm groups that are ecologically char- 
acteristic of TRF can provide important insights into the 
historical construction of the biome [5]. In this respect, 
the pantropical palm family (Arecaceae/Palmae) presents 
an ideal study group. First, palms are among the most 
important and characteristic components of TRF ecosys- 
tems worldwide in terms of species diversity (approxi- 
mately 2,400 species), abundance of individuals and 
impact on the environment [28-30]. Based on the excel- 
lent taxonomic knowledge for this family [28], we calcu- 
lated that over 90% of its species diversity is restricted to 
TRF (Figure 2). Water and energy-related variables are 




Figure 2 Example of an understory lowland tropical rain forest 
in the Parque National do Amazonia (near Itaituba, Para state, 
Brazil) dominated by palms. Foreground Bactris acanthocarpa var. 
exscapa, upper right corner Attalea sp., middle left Euterpe precatoria, 
background: Astrocaryum gynacanthum. Photo: TLPC. 
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strong determinants of palm diversity [31,32] and funda- 
mental anatomical constraints inhibit palms from colo- 
nising cold environments [33,34]. Second, the known 
history of palms extends far back into the Cretaceous 
although the details of the spatiotemporal origin of the 
family remain controversial [28]. Direct evidence from 
unambiguous fossils associated with palms suggest that 
the family was already present during the Turonian (89- 
93.5 Ma, [35-37]) while more doubtful fossils have been 
recorded since the Aptian (112 Ma, [28]). More recently, 
several molecular clock estimates based on angiosperm 
wide phylogenies suggested a stem age for the family ran- 
ging from 91 to 120 Ma [38-41]. These studies were 
based on a very limited sampling within the family and 
thus do not provide reliable approximations for the 
crown node age and early diversification history. To date, 
most estimates of palm ages have focused on subfamily 
[42,43] or tribal levels [44-46]. 

Here we investigate the origin and diversification of 
palms in space and time using the first complete generic- 
level sampling for any important TRF-restricted plant 
family [47] . We estimated speciation events under a Baye- 
sian framework using a relaxed molecular clock approach 
(BEAST, [48]), while the spatial origin of the family was 
inferred under a maximum likelihood method that imple- 
ments the dispersal-extinction-cladogenesis model [49,50]. 
Finally, we use palms as a model to explore the evolution 
of TRF biodiversity by testing which of the three TRF 
diversification hypotheses outlined above corresponds to 
the diversification history of the family. 

Results and discussion 

Evolutionary origin of palms 

The fossil-calibrated molecular dating of a complete 
genus-level palm supertree [47] provides for the first 
time minimum age estimations for all major groups of 
the palm family (Table 1). The resulting chronogram 
(Figure 3b) suggests that the diversification of extant 
lineages of palms started in the mid-Cretaceous period 
at the Albian-Cenomanian boundary (crown node: 100 



Ma, 95% highest posterior density (HPD) 108-92 Ma). 
The Cretaceous represents an important period for 
plant evolution as it witnessed the rise and diversifica- 
tion of flowering plants [51,52]. The fossil record indi- 
cates that from the Albian to the middle Cenomanian 
angiosperms diversified extensively, becoming more 
abundant relative to other plants and establishing them- 
selves as a major part of paleofloras by the end of that 
period [53]. Our maximum likelihood analysis of geo- 
graphic range evolution indicates that the most likely 
distribution of the most recent common ancestor of 
palms was centred on present day Central/North Amer- 
ica and Eurasia, which corresponds to the Laurasian 
landmass at that time (Figure 3a). Notably, the oldest 
reliable palm fossils (Turonian to Campanian) have all 
been discovered in Europe and North America [28]. A 
Laurasian origin for palms was previously suggested by 
Uhl and Dransfield [54] based on the prevalence of 
putatively primitive lineages in the northern hemisphere. 
From this ancestral area several subsequent dispersal 
events are inferred into the equatorial regions of South 
America, Africa and South East Asia, the present day 
distribution of palms. Finally, our ancestral biome analy- 
sis suggests that the earliest palm lineages were 
restricted to TRF (P = 0.984; Figure 3b and Table 1), a 
result that was further supported by a test of phyloge- 
netic signal of the biome. In fact, adaptation to non- 
TRF biomes did not arise until the Paleocene within the 
fan palm subfamily Coryphoideae. Thus, our results sup- 
port the notion that palms originated in a TRF-like 
biome and started to diversify during the mid-Cretac- 
eous period in Laurasia. 

Origin of tropical rain forests 

The presence of TRF during the mid-Cretaceous period 
is controversial because pre-Cenozoic fossils associated 
with TRF are notoriously sparse [8,11,55,56] in contrast 
to the relatively well documented fossil evidence of TRF 
during the early Cenozoic from equatorial to relatively 
high palaeolatitudes [11-14,20,57]. The earliest fossil 



Table 1 Mean estimated ages, 95% confidence intervals and ancestral ecologies for the family and subfamilies 



Clade 


Age in Ma 


95% HPD 


Proportional likelihoods of ancestral 
ecology of branch leading to node 3 


Arecaceae, crown 


100.1 


92-108.7 


0 = 0.984/1 = 0.0/2 = 0.0 


Calamoideae, crown 


80.2 


70.3-90.3 


0 = 0.999/1 = 0.0/2 = 0.0 


Nypoideae, stem 


93.5 


87.5-100.6 


0 = 0.982/1 = 0.0/2 = 0.0 


Coryphoideae, crown 


66.0 


51.35-80 


0 = 0.716/1 = 0.28/2 = 0.0 


Ceroxyloideae, crown 


52.1 


30-74.2 


0 = 0.979/1 = 0.0/2 = 0.0 


Arecoideae, crown 


73.6 


66.1-81.3 


0 = 0.999/1 = 0.0/2 = 0.0 



Bold entries indicate values significantly different from other states. 

a 0 = Tropical rain forest restricted; 1 = mangrove restricted; 2 = non-rain forest restricted. 

HPD = highest posterior density; Ma = millions of years. 
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age in millions of years 

Figure 3 Area, tempo and mode of palm diversification, (a) Paleomap representing the distribution of landmasses in the mid-Cretaceous 
period, dark grey upland land, light grey lowland (100 million years (Ma), adapted from Beerling and Woodward [60]). Laurasia, which is the 
most likely ancestral area reconstructed for the crown node of palms, is highlighted, (b) Chronogram showing the three different biomes 
assigned to each genus. Red: tropical rain forest; green: mangrove; blue: not tropical rain forest; grey: ambiguous. Yellow circles indicate fossil 
calibration points. The vertical black lines highlight the five subfamilies of palms with an illustration (drawings by Marion Ruff Sheehan, LH. Bailey 
Hortorium, Cornell University, except top one (Arecoideae), which is reproduced with permission from Springer from Kahn and de Granville [30]. 
(c) Semilogarithmic mean lineage-through-time (LTT) plot averaged over 1,000 posterior trees from the Bayesian analysis (left axis, triangles) and 
percentage of missing taxa as a function of time (right axis, grey line). Short dashed line = upper 95% confidence interval; long dashed line = 
lower 95% confidence interval; filled square = extant number of palms species. Vertical black line indicates threshold up to which the LTT plot is 
considered reliable even under incomplete taxon sampling. Palm fossil indicates time of earliest known unequivocal fossil for the family 
{Sabalites fossil leaf image reproduced by permission of the Board of Trustees, National Museums Liverpool, Liverpool, UK). 
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flora interpreted as characteristic of TRF was found in 
the early Cenomanian from several formations in North 
America, for example, the Dakota formation of Kansas 
(99 Ma, [9,53,58]), which is consistent with our results. 
However, these conclusions, which were based on the 
physiognomy of leaf characters, such as shape and size, 
that are generally associated with megathermal vegeta- 
tion [20,53], have been questioned by some authors (see 
personal communication from Johnson in Morley [20]). 
Several studies have suggested that during the Cretac- 
eous plant biodiversity was highest at mid paleolatitudes 
where the climate was more favourable while equatorial 
latitudes were exposed to a drier and hotter climate 
unlikely to have supported TRFs [20,59]. In addition, 
simulations of major vegetation distributions during the 
mid-Cretaceous period (100 Ma) indicate that the pre- 
sence of TRF in the Cenomanian of North America and 
Eurasia as well as other parts of the world is plausible 
[14,60]. Finally, indirect evidence is provided by a diver- 
sification study of Malpighiales, a large plant order 
mainly restricted to TRFs. Using a molecular clock 
approach, it was estimated that the origin of this order 
dates to around 114 Ma with subsequent diversification 
during the Cenomanian [10], implying the presence of 
TRF at that time. 

Even though molecular dating methods have been cri- 
ticised in relation to interpretations of TRF origins [12] 
and are not assumption free, such approaches have 
played important roles in understanding the construc- 
tion of other species rich biomes (for example, the Cape 
flora [61,62]) and provide important insights when the 
fossil record is sparse or incomplete. Molecular dating 
of palms, one of the most characteristic TRF plant 
families, provides additional evidence that modern TRF 
might have already been in place 100 Ma, significantly 
earlier than suggested by unequivocal fossil evidence 
[13,57]. It is most likely that formation of the TRF 
biome was a gradual process and thus the precise time 
at which modern TRFs can be recognised may be 
impossible to pinpoint. However, our results and other 
evidence discussed above [9,10,53] imply that the assem- 
bly of the TRF biome had already started during the 
mid-Cretaceous period and was not just a strictly Ceno- 
zoic process. It is puzzling that, to our knowledge, no 
macrofossil of palms has been recovered in Cenomanian 
deposits of North America [9]. However, our results 
suggest that palms were just starting to diversify at this 
time and may not have been widespread, thus reducing 
the frequency of fossilization and probability of later dis- 
covery. Interestingly, studies of other species-rich TRF 
plant families yield timings for the earliest extant lineage 
diversification events (that is, crown node estimates) 
that largely post-date palms, for example Leguminosae 
(59 Ma, [63]), Annonaceae (89 Ma, [64]) and Rubiaceae 



(86 Ma, [65]). This would imply that palms represent 
one of the first extant plant families to have diversified 
within TRF since its origin. Thus, palms not only play a 
major role in present day TRF [28-30], but also appear 
to have been a key component in the assembly and 
diversity of this biome since the earliest stages of its 
evolution. 

Early diversification of tropical rain forests 

To depict global diversification at the family level and 
test the hypotheses of TRF diversification, we generated 
a semilogarithmic lineage-through-time (LTT) plot as 
well as 95% confidence intervals based on 1,000 ran- 
domly selected posterior trees from the BEAST output 
(Figure 3c). LTT plots are widely used to characterise 
the diversification of clades as a function of time 
[66-68], but they are sensitive to incomplete taxon sam- 
pling [69,70], as is the case here (183 species sampled 
out of circa 2,400, approximately 7.6%, see Additional 
file 1). However, our sampling is phylogenetically repre- 
sentative and non-random as we included 100% of all 
described palm genera. Such a sampling strategy has the 
advantage of representing all the deeper nodes of the 
phylogeny and can provide a good estimation of diversi- 
fication history up to a point in time after which under- 
sampling at shallower nodes biases the inference [71,72]. 
In order to restrict the interpretation of our LTT plot to 
the more accurately estimated parts we used a novel 
approach to identify a threshold after which the poten- 
tial impact of incomplete taxon sampling becomes too 
important to permit accurate analyses of diversification 
rates (see Methods). Based on the chronogram, the per- 
centage of total missing palm species increased sharply 
from 12% to 28% after circa 24 Ma (Figure 3c). Interpre- 
tation of the LTT plot was thus restricted to the period 
prior to that point in time (from 100 Ma to 24 Ma) and 
all nodes occurring after the 24 Ma threshold were 
excluded from subsequent analyses. 

Between 100 Ma and 24 Ma, the LTT plot forms an 
almost straight line, which suggests that the palm family 
underwent a constant rate of diversification without any 
major shifts or radiations [67,68]. This is also statisti- 
cally supported by the better fit of the pure-birth model 
of diversification (constant diversification rates with no 
extinction) on the LTT plot than any other model tested 
(AAIC RC = -1.421 between the pure birth model and the 
second best fitting model (the density-dependent or 
'DDX' model); AIC = Akaike Information Criterion). 
Moreover, the AAIC RC was not significantly different 
under the null hypothesis of diversification rate consis- 
tency when calculated from 10,000 phylogenies simu- 
lated under the pure-birth process (P value = 0.739). 
Finally, the y statistic [70] also supported a constant rate 
hypothesis as it was not significantly different from zero 
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(y = -0.318, P value = 0.379). These results lend support 
to the museum model (Figure 1) in which diversification 
rates remain constant and extinction rates are low 
[16,22]. This hypothesis was thought to be the conse- 
quence of the old and ecologically stable conditions of 
TRF over millions of years. Even though such views 
have now been replaced by the notion of greater dyna- 
mism in TRF evolution [14,17,73], the biome itself has 
never completely disappeared [14,20] and has persisted 
in refuges during unfavourable climatic times. The exis- 
tence of such refuges may have allowed comparatively 
old lineages to persist and contribute to present day 
species diversity. For example, lineages of another 
diverse TRF plant family, Annonaceae, were shown to 
have persisted in possible TRF refuges of East and West 
Africa for over 30 million years [74] even during clima- 
tically unfavourable times (for example, the global cool- 
ing of the Eocene/Oligocene boundary). We suggest 
here that TRF refugia may have played a similar role 
throughout the history of the palm family and, as a 
result, global palm species diversity is at least partly the 
result of a gradual accumulation of ancestral lineages 
through time, and cannot be attributed to ancient or 
recent speciation bursts alone. 

At finer time scales, palm species diversity and diversi- 
fication rates most likely fluctuated with extinction rates 
possibly increasing and decreasing at specific time 
frames in the past, perhaps in relation to climatic and 
geological changes. For example, studies of the palyno- 
flora through the Late Paleocene-Eocene Thermal Maxi- 
mum in Colombia indicate an increase in palm 
morphospecies after this time (56.3 Ma, [8]). Thus, it is 
probable that different clades within palms underwent 
alternative diversification scenarios leading to a hetero- 
geneous pattern of diversification among the lineages 
within the family. However, our study implies that over 
a larger time scale these changes did not influence the 
overall global pattern of diversification in palms, at least 
until 24 Ma. This is also apparent during the major 
extinction episode at the Cretaceous/Paleogene bound- 
ary (K/Pg; [75]), which had no statistically significant 
effect on diversification rates (Figure 3c). 

To date, relatively few family-level studies have pro- 
vided evidence for the museum model of tropical plant 
diversification (Annonaceae [64], liverwort family Lejeu- 
naceae [76]). This pattern contrasts with the study of 
Malpighiales evolution, which indicated that all major 
lineages originated within a short timeframe suggesting 
an early rapid speciation of the order, although no 
detailed diversification analysis was undertaken for this 
group [10]. Interestingly, meta-analyses based on a large 
number of species-level dated molecular phylogenies of 
a range of plants and animals have also underlined the 
importance of the constant rate diversification model 



[77,78]. Indeed a large number of phylogenies fitted the 
simplest model of diversification. For example, Morion 
et al. [78], using a novel coalescent-based approach, 
found that 87 out of 289 phylogenies studied (30%) bet- 
ter fitted the Yule process (time constant rates with no 
extinction) than alternative models. Even though these 
results were obtained from a wide range of organisms 
occurring in different ecosystems, it nevertheless under- 
lines the importance of such a simple process for 
explaining present day diversity. 

Given the threshold of 24 Ma imposed on our LTT 
plot (Figure 3c), there is little indication about recent 
(Neogene) diversification patterns. However, it is clear 
that in order to attain present day species diversity 
(Figure 3c), rates must have increased, which suggests a 
shift in diversification that occurred after 24 Ma. This 
could have been achieved either by accelerating diversi- 
fication rates, consistent with the 'recent cradle' model, 
or simply by higher overall constant rates. Recent high 
speciation rates within other TRF plant genera have 
been documented [13-15] and it seems likely that rapid 
speciation occurred within some species-rich palm gen- 
era given the very young stem node age estimations we 
inferred for them (for example, Pinanga, approximately 
130 species, stem node approximately 12 Ma; Dypsis, 
approximately 140 species, stem node approximately 13 
Ma). In fact, the 'museum' and 'recent cradle' hypoth- 
eses are not mutually exclusive; both mechanisms could 
be at work, possibly within different palm lineages and 
different time frames. For example, both the 'ancient' 
and 'recent' cradle models of diversification have 
been identified within the TRF-restricted leaf beetle 
family [79]. 

Conclusions 

Our results from one of the most important TRF plant 
families suggests that present day TRF biodiversity can 
at least in part be explained by a steady accumulation 
of lineages dating back to the mid-Cretaceous period 
and is not just the result of rapid radiations, either 
recent [25] or ancient [10]. The analysis of additional 
family-level diversification patterns of other TRF 
restricted plant groups will undoubtedly shed more 
light into the evolutionary forces that have led to the 
immense diversity of species found within modern TRF 
today. 

Methods 

Taxon sampling 

This study builds upon the complete generic-level 
supertree analyses of palms by Baker et al. [47], the 
most extensive phylogenetic study of the family pub- 
lished to date. Here, the sampling is updated to be con- 
sistent with the latest family- wide monograph [28,80], 
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in which 183 genera are accepted (see Additional file 1 
for the list of genera used). This was performed by 
repeating the supertree analyses of Baker et al. [47] 
with the addition of published plastid DNA sequence 
data for the recently described genus Tahina [81] (see 
below). Three commelinid monocot outgroup taxa (Cos- 
tus, Dasypogon, Zed) were selected from the sampling 
of Baker et al. [47]. 

Fossil calibration 

Palms have a rich fossil record dating from the Late 
Cretaceous onwards. Although the record is unusually 
rich among angiosperms, only a small fraction of palm 
fossils can be identified to specific taxonomic groups 
with confidence. Drawing on recent surveys of the palm 
fossil record [28,82,83], we selected the most reliable 
fossils (Table 2), judged on the basis of the credibility of 
their purported taxonomic affinities and reported ages. 
Nevertheless, none of these fossils is sufficiently infor- 
mative to justify allocation to crown nodes [84]; they are 
thus applied conservatively to stem nodes throughout. 
Where authors provide a range of age estimates, we 
have used the most recent date. Where a geological 
time period alone is specified, we have used the date of 
the upper end of that period [85]. 

The four selected fossils are widely distributed across 
the family and are located in three out of the five subfa- 
milies. The earliest fossils that can be assigned unequi- 
vocally to a taxonomic group within palms are Late 
Cretaceous records of palmate leaves, the earliest of 
which is Sabalites carolinensis from the late Coniacian 
of South Carolina [86] . Although an affinity with Sabal 
is implied by the genus name, the fossil could be linked 
with many coryphoid groups and its age is therefore 
used conservatively here as a calibration point for the 
stem node of subfamily Coryphoideae as a whole with 
an age of 85.8 Ma. Hyphaene kapelmanii, a fossil discov- 
ered at a late Oligocene site in Ethiopia [87], provides a 
further calibration point within the Coryphoideae. This 
fossil consists of a petiole fragment with a close resem- 
blance to the modern genus Hyphaene due to the mor- 
phology of its wide, recurved marginal spines. We use 



Table 2 Names of fossils used to calibrate the tree, with 
the respective exponential prior parameters used 



Fossil name 


Hard lower 
bound (Ma) 


Soft upper 
bound 95% (Ma) 


Exponential mean 
(uncertainty) 


Sabalites 
carolinensis 


85.8 


88.8 


1 


Mauritiidites 


65 


69.49 


1.5 


Attaleinae 


54.8 


60.79 


2 


Hyphaene 
kapelmanii 


27 


28.5 


0.5 



Ma = millions of years. 



this fossil as a constraint for the stem node of subtribe 
Hyphaeninae with an age of 27 Ma. 

In subfamily Calamoideae, the unique structure of the 
pollen of subtribe Mauritiinae corresponds closely to 
fossil pollen in the genus Mauritiidites, specifically the 
clavate monosulcate grains with each spine inset and a 
swollen foot layer below. Mauritiidites has been 
recorded as early as the Maastrichtian of Africa [88] 
with numerous records soon after in the Palaeocene 
onwards of South America [89]. We use it here as a 
calibration for the stem node of the Mauritiinae with an 
age of 65 Ma. 

Fossil records of the coconut tribe Cocoseae, particularly 
of fossilised endocarps, are numerous [28] . Until recently, 
well documented records appeared from the Middle 
Eocene onwards (for example, [90,91]), but new research 
in the middle to late Palaeocene of Colombia has revealed 
compression fossils of large fruits that closely resemble the 
modern coconut, Cocos nucifera, both in size and surface 
morphology [92]. In the absence of further substantiating 
evidence, we allocate this fossil to the stem node of the 
Attaleinae, the subtribe of tribe Cocoseae to which Cocos 
belongs, with an age of 54.8 Ma. 

Finally, a number of other reliable fossils could not be 
used because they are made redundant by older fossils 
assigned to more distal nodes. Nypa is most notable 
here, given its outstanding macrofossil and microfossil 
records dating back to the Maastrichtian [28,93]. Also 
significant are the distinctive diaperturate fossil pollen 
grains, usually referred to the form genus Dicolpopollis, 
which can be assigned with confidence to tribe Cala- 
meae of the Calamoideae [82,83,94-96]. The earliest 
records of this fossil palynomorph are from the Maas- 
trichtian and its boundary with the early Palaeocene of 
Somalia and Borneo [88,96]. 

For each fossil we applied an exponential prior, the 
parameters of which are given in Table 2. Finally, the 
stem node of palms was constrained by a uniform prior 
ranging from 110 to 120 Ma. This corresponds to the 
earliest monocot fossil [97]. By doing this we imply that 
the stem of palms cannot be older than the oldest 
monocot fossil. 

Molecular dating 

Molecular dating was carried out using BEAST 1.5.3 
[48,98]. For this analysis the 'most congruent supertree', 
based on the method and data of Baker et al. [47], was 
used as a topological constraint. This topology was 
based on extensive data sampling, including DNA 
sequence data, restriction fragment length polymorph- 
isms (RFLP) and morphology, and represents the best 
family-wide estimation of phylogenetic relationships 
between palm genera to date. To update the taxonomic 
sampling of the supertree, we repeated the supertree 
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analysis of Baker et al. [47]. This study used a matrix 
representation with parsimony (MRP) analysis based on 
input trees generated from individual partitions and 
combinations of partitions with matrix elements 
weighted in proportion to bootstrap values of corre- 
sponding input tree clades. The strict consensus tree of 
this analysis was highly resolved with minor ambiguity 
in tribe Trachycarpeae and parts of tribe Areceae only. 
One most-parsimonious tree was selected at random 
and used as a constraint. This tree was pruned to 
include only the 183 genera accepted by Dransfield et 
al. [28]. All DNA sequence datasets utilised by Baker et 
al. [47] were included in our molecular dating analysis 
(plastid DNA regions: atpB, matK, ndhF, rbcL, rpsl6 
intron, trnD-trnT, trnL-trnF, trnQ-rpsl6; nuclear DNA 
regions: 18S, ITS, ms, prk, rpb2). The completeness of 
taxonomic sampling for each of these regions varies 
from 12% to 100%, with an average of 48%. Moreover, 
sampling for the chloroplast markers was much more 
complete than for the nuclear markers. Morphological 
and RFLP datasets used by Baker et al. [47] were 
excluded. The supertree topology was used as constraint 
by deleting in the XLM BEAST input file the following 
commands: subtreeSlide; narrowExchange; wideEx- 
change; wilsonBalding. Each marker was individually 
partitioned in BEAUTi 1.5.3 http://beast.bio.ed.ac.uk/ 
BEAUti and assigned the General Time Reversible 
model (GTR) with gamma-distributed rate variation (G). 
model of sequence evolution. Prior to our full analysis, 
we investigated the effect of missing data on the estima- 
tion of ages. We undertook a preliminary analysis on 2 
datasets: 1 with all the 13 markers (missing data pre- 
sent) and 1 where the 5 nuclear markers were removed 
(missing data of 10%). A regression analysis between the 
ages obtained for all nodes was highly significantly posi- 
tive (R = 0.88; t test: P < 0.001) indicating that missing 
data in our dataset are likely to have little influence on 
age estimations in molecular dating. We then undertook 
the full-scale analysis on the full 13-marker dataset. 

In total, 8 individual analyses were carried out, 5 with 
20 million generations and 2 with 30 million genera- 
tions, resulting in a total of 160 million generations, and 
sampling every 1,000th generation. Individual analyses 
were performed in order to test for convergence of the 
results. Tracer 1.4 [99] was used to check for conver- 
gence of the model likelihood and parameters between 
each run. Results were considered reliable once the 
effective sampling sizes (ESS) of all parameters exceeded 
200. The resulting independent log and tree files were 
then combined using LogCombiner discarding 10% of 
generations as burn in per independent run. 

Finally, deviation from a strict molecular clock was 
tested by running the analysis a second time with the 
strict clock enforced. We used the Bayes Factor as 



implemented in Tracer 1.4 [99] to select the best-fitting 
model under the smoothed marginal likelihood estimate 
and with 100 bootstrap replicates [100]. This test 
strongly supported the data as being non-clock like (In 
BF = 1,158.2 ± 3.2 in favour of relaxed clock hypoth- 
esis), and thus the results under the relaxed clock are 
presented here. 

Diversification analyses 

To explore diversification rates in the family, we gener- 
ated a semilogarithmic LTT plot. The mean LTT plot as 
well as the 95% confidence intervals was generated from 
a random selection of 1,000 posterior trees resulting 
from the BEAST analysis. LTT plots are sensitive to 
incomplete taxon sampling [70]. However, the full gen- 
eric sampling of palm genera (100% of genera included) 
means that sampling of extant lineages is complete in 
the older parts of the phylogeny, becoming progressively 
more incomplete towards the present [71,101]. In order 
to avoid misinterpretation of the plot we restrict our 
analyses to the accurately estimated part of the LTT, 
which will be least influenced by the missing taxa. We 
do this by finding the point at which incomplete taxon 
sampling will likely begin to have a significant effect on 
the LTT plot. Under the assumption that all genera are 
monophyletic, speciation events within each genus will 
always be found after the stem node of that genus. 
Given the known age of each stem node for each genus, 
we calculated a cumulative total for the number of miss- 
ing species as a function of time. By doing this we gen- 
erated a time-dependent curve representing the increase 
of missing taxa from the origin of the family until the 
present that quantified the amount of uncertainty each 
part of the LTT plot contains (Figure 3c). In this study, 
the proportion of missing taxa was less than 12% from 
100 Ma to 24 Ma, at which time a dramatic increase 
occurred with missing species rising to 24% and on to > 
92% at the present time. This point represents the stem 
age of the genus Calamus, the most species-rich genus 
in palms [28]. All nodes occurring after this 24 Ma 
threshold were deleted from subsequent diversification 
analyses. 

Two different approaches were used to test for signifi- 
cant changes in diversification rates. First, we used a 
maximum likelihood method for fitting alternative diver- 
sification models to the LTT plot [102] using the R 
package LASER 2.2 [103]. The test statistic for diversifi- 
cation rate constancy AAIC RC is calculated as AAIC RC = 
AIC R c - AIC Rv , where AIC Rc is the AIC score for the 
best fitting rate-constant diversification model, and 
AICr v is the AIC for the best fitting variable-rate diver- 
sification model. A negative value for AAIC RC indicates 
that the data is best approximated by a rate-constancy 
model. We fitted five different diversification models: (1) 



Couvreur et al. BMC Biology 201 1, 9:44 
http://www.biomedcentral.eom/1 741 -7007/9/44 



Page 9 of 12 



the constant-rate birth model (the Yule process; [104]) 
with the speciation rate (k) being constant and the 
extinction (u) set to zero; (2) the constant-rate birth- 
death model with two parameters, speciation (k) and 
extinction (u); (3) a pure birth rate-variable model 
where the speciation rate XI shifts to rate X2 at time ts, 
with three parameters (kl, X2, ts); (4) an exponential 
density-dependent speciation rate 'DDX' model; and (5) 
a logistic density-dependent speciation rate 'DDL' 
model. The significance of the observed AAIC RC was 
evaluated by simulating 10,000 trees under a pure birth 
constant diversification rate. 

Second, we calculated the y statistic of Pybus and Har- 
vey [70], which provides a summary of the distribution 
of nodes in the phylogeny: if the internal nodes are clo- 
ser to the root then y < 0; if they are closer to the tips 
then y > 0; if the nodes are equally spread out then y = 
0. The observed y statistic was compared with the distri- 
bution of the y statistic of 1,000 simulated phylogenies 
under a pure-birth model using the LASER package 2.2. 

Ancestral areas 

A presence-absence matrix was built representing the 
global distribution of palm genera (see Additional file 1 
for the original data used to perform this analysis). We 
defined seven, non-overlapping major palm areas that 
reflect the distribution and endemism of genera as well 
as broad scale geological units, as follows: (A) South 
America, (B) North America (including Central Ameri- 
can and the Caribbean), (C) Africa (including Arabia); 
(D) Indian Ocean (Madagascar, Mascarenes, Comoros 
and Seychelles), (E) India (including Sri Lanka), (F) Eur- 
asia (including west Malesia to the west of Wallace's 
Line) and (G) Pacific Ocean (including east Malesia to 
the east of Wallace's Line, Australia and the Pacific 
Islands). Each genus was assigned to one or more of the 
major palm areas based on its known current distribu- 
tion [28,105]. 

Ancestral areas (AA) were reconstructed using a maxi- 
mum likelihood method under the dispersal-extinction- 
cladogenesis model [49,50] as implemented in the soft- 
ware Lagrange build 20091004 [50]. We tested our 
ancestral area reconstruction under two different bio- 
geographic models (see Additional file 1 for the para- 
meters used to perform these analyses). The first model 
(M 0 ) was unconstrained and we assigned an equal prob- 
ability (P = 1.00) of dispersal between all areas during 
the whole time period considered. This model assumes 
that spatial relationships among areas have no effect on 
biogeographical patterns. For the second model (Mi), 
we applied a more complex biogeographic scenario 
incorporating prior information on range evolution as 
well as dispersal probabilities between areas given dis- 
crete periods of time. This model was based on past 



climatic data, tectonic history and presence/absence of 
postulated land bridges [14,19,20,106-108]. Five time 
frames were delimited and dispersal probabilities were 
assigned between all adjacent areas (see Additional file 
1). Dispersal probabilities were set as following: low or 
no dispersal = 0.01; low dispersal = 0.25; medium dis- 
persal = 0.5; high dispersal = 0.75; areas adjacent or very 
close = 1. Scripts in the programming language Python 
http://www.python.org/ were generated using the online 
helper http://www.reelab.net/lagrange. Because of the 
large surface of each area (continent-level areas), the 
maximum number of ancestral areas was limited to two. 

The genus Cocos presented a special problem because 
it is widely distributed across all areas. Such highly poly- 
morphic states generally inject a high level of ambiguity 
into the analyses, as was the case in preliminary analyses 
here. Following the recommendations of Ronquist [109], 
we allocated a putative ancestral area coding to Cocos 
(area A) based on the findings of Meerow et al. [110]. 

Ancestral biome 

The probable ancestral biome at the crown node of 
palms was reconstructed under a maximum likelihood 
method using the Markov k-state 1 parameter model 
(MkjJ model of character evolution implemented in 
Mesquite version 2.74 [10,111]. We assigned genera into 
three different biome state categories following Olson 
et al. [112]. State 0: predominantly 'tropical and subtro- 
pical moist broadleaf forests' biome (that is, TRF); state 
1: 'mangrove'; state 2: a general category that contains 
genera not belonging to any of the two first categories 
(that is, not TRF-restricted). State 2 is broad in its defi- 
nition encompassing all other biomes for palms. When 
a genus occurred in both TRF and non-TRF biomes it 
was coded as ambiguous (see Additional file 1 for the 
original data used to perform this analysis). The man- 
grove category was included to take into account the 
ecology of Nypa. The phylogenetic signal of the biome 
character was tested by randomising the tips of the phy- 
logeny 1,000 times in Mesquite [111] in order to create 
a null distribution of the number of steps under the 
maximum parsimony criterion and the 99% confidence 
intervals. This null distribution was compared to the 
observed number of steps necessary to explain the 
occurrence of each character on the phylogeny. In this 
case, the observed value fell outside the 99% confidence 
interval confirming that the biome category is phylogen- 
etically conserved. 

Additional material 



Additional file 1: Additional tables. Table SI: Genera sampled with 
total number of species per genus, biome coding used in Mesquite and 
area coding used in the Lagrange analysis. Table S2: Alternative dispersal 
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models between areas used in Lagrange. This file presents the names of 
all officially recognised palm genera, with the coding for present day 
biome and area, as well as the details of the two alternative 
biogeographical models used in the analysis. 



Acknowledgements 

We are grateful to Jim Doyle, Roy Erkens, Alan Graham, Madeline Hariey, 
David Mabberley, Jim Miller, Bob Morley, Vincent Savolainen, Selena Smith 
and Gary Upchurch for their advice during the preparation of this 
manuscript. Alexander Antonelli and one anonymous reviewer are gratefully 
acknowledged for their comments and suggestions. Gregory Plunkett and 
Ken Carol kindly provided logistical support for our analyses. Tylor Durand is 
thanked for his help in preparing the photo of Figure 2. This work was 
supported by a grant from the UK Natural Environment Research Council to 
WJB. 

Author details 

'The New York Botanical Garden, 200th Street and Kazimiroff Boulevard, 
Bronx, NY 10458-5126, USA. 2 Royal Botanic Gardens, Kew, Richmond, Surrey, 
TW9 3AB, UK. 3 lnstitut de Recherche pour le Developpement (IRD), UMR 
DIA-DE, DYNADIV researche group, 91 1, avenue Agropolis, BP 64501, F-34394 
Montpellier cedex 5, France. 

Authors' contributions 

TLPC and WJB designed the study; TLPC, FF and WJB performed the 
research, analysed the data and wrote the paper. All authors read and 
approved the final manuscript. 

Received: 8 April 201 1 Accepted: 16 June 201 1 Published: 16 June 201 1 
References 

1 . Gentry A: Tropical forest biodiversity: distributional patterns and their 
conservational significance. Oikos 1992, 63:19-28. 

2. Whitmore TC: An introduction to tropical rain forests Oxford, UK: Clarendon 
Press; 1998. 

3. Hill JL, Hill RA: Why are tropical rain forests so species rich? Classifying, 
reviewing and evaluating theories. Prog Phys Geog 2001, 25:326-354. 

4. Asner GP, Rudel TK, Aide TM, Defries R, Emerson R: A contemporary 
assessment of change in humid tropical forests. Conserv Biol 2009, 
23:1386-1395. 

5. Pennington RT, Richardson JE, Lavin M: Insights into the historical 
construction of species-rich biomes from dated plant phylogenies, 
neutral ecological theory and phylogenetic community structure. 

New Phytol 2006, 172:605-616. 

6. Moritz C, Patton JL, Schneider CJ, Smith TB: Diversification of rainforest faunas: 
an integrated molecular approach. Annu Rev Ecol Syst 2000, 31:533-563. 

7. Jaramillo C, Rueda MJ, Mora G: Cenozoic plant diversity in the neotropics. 
Science 2006, 311:1893-1896. 

8. Jaramillo C, Hoorn MC, Silva S, Leite F, Herrera F, Quiroz L, Dino R, 
Antonioli L: The origin of the modern Amazon rainforest: implications 
from the palynological and paleobotanical record. In Amazonia, 
Landscape and Species Evolution. Edited by: Hoorn MC, Wesselingh FP. 
Oxford, UK: Blackwells; 2010:317-334. 

9. Upchurch GR, Wolf JA: Mid-Cretaceous to early tertiary vegetation and 
climate: evidence from fossil leaves and woods. In The Origin of 
Angiosperms and Their Biological Conseguences. Edited by: Friis EM, Chaloner 
WG, Crane PR. Cambridge, UK Cambridge University Press; 1987:75-105. 

10. Davis CC, Webb CO, Wurdack KJ, Jaramillo CA, Donoghue MJ: Explosive 
radiation of malpighiales supports a mid-Cretaceous origin of modern 
tropical rain forests. Am Nat 2005, 165:E36-E65. 

1 1. Johnson KR, Ellis B: A tropical rainforest in Colorado 1.4 million years 
after the Cretaceous-Tertiary boundary. Science 2002, 296:2379-2383. 

12. Wing SL, Herrera F, Jaramillo CA, Gomez-Navarro C, Wilf P, Labandeira CC: 
Late Paleocene fossils from the Cerrejon formation, Colombia, are the 
earliest record of neotropical rainforest. Proc Natl Acad Sci USA 2009, 
106:18627-18632. 



13. Jacobs BF: Palaeobotanical studies from tropical Africa: relevance to the 
evolution of forest, woodland and savannah biomes. Philos Trans ft Soc 
Lond, B 2004, 359:1573-1583. 

14. Morley RJ: Origin and Evolution of Tropical Rain Forests New York, USA: John 
Wiley & Sons; 2000. 

15. Raven HP, Axelrod Dl: Angiosperm biodiversity and past continental 
movements. Ann Missouri Bot Gard 1974, 61:539-673. 

16. Wallace AR: Tropical Nature, and Other Essays London, UK Macmillian; 1878. 

17. Richards PW: The Tropical Rain Forest: an Ecological Study. 2 edition. 
Cambridge, UK Cambridge University Press; 1996. 

18. Wolfe JA: Some aspects of plant geography of northern hemisphere 
during Late Cretaceous and Tertiary. Ann Missouri Bot Gard 1975, 
62:264-279. 

19. Zachos J, Pagani M, Sloan L, Thomas E, Billups K: Trends, rhythms, and 
aberrations in global climate 65 Ma to present. Science 2001, 292:686-693. 

20. Morley R: Cretaceous and Tertiary climate change and the past 
distribution of megathermal rainforests. In Tropical Rainforest Responses to 
Climatic Changes. Edited by: Bush MB, Flenley J. Chichester, UK: Praxis 
Publishing; 2007:1-31. 

21 . Wilf P, Cuneo NR, Johnson KR, Hicks JF, Wing SL, Obradovich JD: High 
plant diversity in Eocene South America: evidence from Patagonia. 
Science 2003, 300:122-125. 

22. Stebbins GL: Flowering Plants: Evolution Above the Species Level Cambridge, 
MA: Harvard University Press; 1974. 

23. Haffer J: Speciation in Amazonian forest birds. Science 1969, 165:131-137. 

24. Gentry A: Neotropical floristic diversity: phytogeographical connections 
between Central and South America, Pleistocene climatic fluctuations, or 
an accident of the Andean orogeny? Ann Missouri Bot Gard 1982, 
69:557-593. 

25. Richardson JE, Pennington RT, Pennington TD, Hollingsworth PM: Rapid 
diversification of a species-rich genus of neotropical rain forest trees. 
Science 2001, 293:2242-2245. 

26. Janssens SB, Knox EB, Huysmans S, Smets EF, Merckx V: Rapid radiation of 
Impatiens (Balsaminaceae) during Pliocene and Pleistocene: result of a 
global climate change. Mol Phylogenet Evol 2009, 52:806-824. 

27. Erkens RHJ, Chatrou LW, Maas JW, van der Niet T, Savolainen V: A rapid 
diversification of rainforest trees (Guatteria; Annonaceae) following 
dispersal from Central into South America. Mol Phylogenet Evol 2007, 
44:399-411. 

28. Dransfield J, Uhl NW, Asmussen CB, Baker JW, Hariey MM, Lewis CE: Genera 
Palmarum: The Evolution and Classification of Palms Kew, UK Kew 
Publishing; 2008. 

29. Phillips O, Miller JS: In Global Patterns of Plant Diversity: Alwyn H. Gentry's 
Forest Transect Data Set. Volume 84. Louis, MO: Missouri Botanical Press; 
2002(St). 

30. Kahn F, de Granville J-J: Palms in Forest Ecosystems of Amazonia Berlin, 
Germany: Springer-Verlag; 1992. 

31. Bjorholm S, Svenning JC, Skov F, Balslev H: Environmental and spatial 
controls of palm (Arecaceae) species richness across the Americas. 
Glob Ecol Biogeogr 2005, 14:423-429. 

32. Kreft H, Sommer JH, Barthlott W: The significance of geographic range 
size for spatial diversity patterns in neotropical palms. Ecography 2006, 
29:21-30. 

33. Tomlinson PB: The uniqueness of palms. Bot J Linn Soc 2006, 151:5-14. 

34. Tomlinson PB: Systematics and ecology of the Palmae. Annu Rev Ecol Syst 
1979, 10:85-107. 

35. Kvacek J, Herman AB: Monocotyledons from the early Campanian 
(Cretaceous) of Grunbach, lower Austria. Rev Paleobot Palyno 2004, 
128:323-353. 

36. Daghlian CP: A review of the fossil record of monocotyledons. Sot Rev 

1981, 47:517-555. 

37. Crie L: Studies on the silicified palms of the Cretaceous formation of 
Anjou [translated from French]. Bull Soc Etudes Sci Angers 1892, 21:97-103. 

38. Bremer K: Early Cretaceous lineages of monocot flowering plants. Proc 
Natl Acad Sci USA 2000, 97:4707-471 1 . 

39. Anderson CL, Janssen T: Monocots. In The TimeTree of Life. Edited by: 
Hedges BS, Kumar S. New York, NY: Oxford University Press; 2009:203-212. 

40. Janssen T, Bremer K: The age of major monocot groups inferred from 
800 + rbcL sequences. Bot J Linn Soc 2004, 146:385-398. 

41 . Wikstrdm N, Savolainen V, Chase MW: Evolution of the angiosperms: 
calibrating the family tree. Proc Natl Acad Sci USA 2001, 268:221 1-2220. 



Couvreur et al. BMC Biology 201 1, 9:44 
http://www.biomedcentral.eom/1 741 -7007/9/44 



Page 11 of 12 



42. Savolainen V, Anstett MC, Lexer C, Hutton I, Clarkson JJ, Norup MV, 
Powell MP, Springate D, Salamin N, Baker WJ: Sympatric speciation in 
palms on an oceanic island. Nature 2006, 441:210-213. 

43. Trenel P, Gustafsson MHG, Baker WJ, Asmussen-Lange CB, Dransfield J, 
Borchsenius F: Mid-Tertiary dispersal, not Gondwanan vicariance explains 
distribution patterns in the wax palm subfamily (Ceroxyloideae: 
Arecaceae). Mol Phylogenet Evol 2007, 45:272-288. 

44. Gunn FB: The phylogeny of the Cocoeae (Arecaceae) with emphasis on 
Cocos nucifera. Ann Missouri Sot Card 2004, 91:505-522. 

45. Cuenca A, Asmussen-Lange CB, Borchsenius F: A dated phylogeny of the 
palm tribe Chamaedoreeae supports Eocene dispersal between Africa, 
North and South America. Mol Phylogenet Evol 2008, 46:760-775. 

46. Roncal J, Borchsenius F, Asmussen-Lange CB, Balslev H: Divergence times 
in the tribe Geonomateae (Arecaceae) coincide with Tertiary geological 
events. In Diversity, Phylogeny and Evolution of Monocotyledons. Edited by: 
Seberg 0, Peterson AT, Barfod A, Davis AM. Aarhus, Denmark: Aarhus 
University Press; 2010:245-265. 

47. Baker WJ, Savolainen V, Asmussen-Lange CB, Chase MW, Dransfield J, 
Forest F, Harley MM, Uhl NW, Wilkinson M: Complete generic-level 
phylogenetic analyses of palms (Arecaceae) with comparisons of 
supertree and supermatrix approaches. Syst Biol 2009, 58:240-256. 

48. Drummond AJ, Ho SYW, Phillips MJ, Rambaut A: Relaxed phylogenetics 
and dating with confidence. P/oS Biol 2006, 4:699-710. 

49. Ree RH, Moore BR, Webb CO, Donoghue MJ: A likelihood framework for 
inferring the evolution of geographic range on phylogenetic trees. 
Evolution 2005, 59:2299-2311. 

50. Ree RH, Smith SA: Maximum likelihood inference of geographic range 
evolution by dispersal, local extinction, and cladogenesis. Sysf Biol 2008, 
57:4-14. 

51. Crane PR, Friis EM, Pedersen KR: The origin and early diversification of 
angiosperms. Nature 1995, 374:27-33. 

52. Friis EM, Pedersen KR, Crane PR: When Earth started blooming: insights 
from the fossil record. Curr Opin Plant Biol 2005, 8:5-12. 

53. Upchurch GR, Wolfe JA: Cretaceous vegetation of the western interior 
and adjacent regions of North America. In Evolution of the Western Interior 
Basin. Volume 39. Edited by: Kaufman EG, Caldwell WGE. St. John's, Canada: 
Geological Association of Canada; 1993:243-281. 

54. Uhl NW, Dransfield J: Genera Palmarum Lawrence, KS, USA: Allen Press; 1 987. 

55. Tiffney BH: Seed size, dispersal syndromes, and the rise of the 
Angiosperms: evidence and hypothesis. Ann Missouri Bot Gard 1984, 
71:551-576. 

56. Wing SL, Boucher LD: Ecological aspects of the Cretaceous flowering 
plant radiation. Ann Rev Earth Planet Sci 1998, 26:379-421. 

57. Burnham RJ, Johnson KR: South American palaeobotany and the origins 

of neotropical rainforests. Philos Trans ft Soc land, B 2004, 359:1595-1610. 

58. Wolfe JA, Upchurch GR: North American nonmarine climates and 
vegetation during the late Cretaceous. Palaeogeogr Palaeoclimatol 
Palaeoecol 1987, 61:33-77. 

59. DiMichele WA, Gastaldo RA: Plant paleoecology in Deep Timel. Ann 
Missouri Bot Gard 2008, 95:144-198. 

60. Beerling DJ, Woodward Fl: Vegetation and the Terrestrial Carbon Cycle: 
Modelling the First 400 Million Years Cambridge, UK: Cambridge University 
Press; 2001. 

61. Linder HP: The radiation of the Cape flora, southern Africa. Biol Rev 2003, 
78:597-638. 

62. Linder HP: Evolution of diversity: the Cape flora. Trends Plant Sci 2005, 
10:536-541. 

63. Lavin M, Herendeen PS, Wojciechowski MF: Evolutionary rates analysis of 
Leguminosae implicates a rapid diversification of lineages during the 
tertiary. Sysf Biol 2005, 54:575-594. 

64. Couvreur TLP, Pirie MD, Richardson JE, Chatrou LW, Saunders RMK, Su YCF, 
Erkens RHJ: Early evolutionary history of the flowering plant family 
Annonaceae: steady diversification and boreotropical geodispersal. 

J Biogeogr 2011, 38:664-680. 

65. Bremer B, Eriksson T: Time tree of Rubiaceae: phylogeny and dating the 
family, subfamilies, and tribes. Int J Plant Sci 2009, 170:766-793. 

66. Ricklefs RE: Estimating diversification rates from phylogenetic 
information. Trends Ecol Evol 2007, 22:601-610. 

67. Nee S, May RM, Harvey PH: The reconstructed evolutionary process. 
Phil Trans R Soc Lond B 1 994, 344:305-31 1 . 



68. Nee S, Mooers AO, Harvey PH: Tempo and mode of evolution revealed 
from molecular phytogenies. Proc Natl Acad Sci USA 1992, 89:8322-8326. 

69. Harvey P, Pagel M: The Comparative Method in Evolutionary Biology Oxford, 
UK: Oxford University Press; 1991. 

70. Pybus OG, Harvey PH: Testing macro-evolutionary models using 
incomplete molecular phytogenies. Proc ft Soc B Biol Sci 2000, 
267:2267-2272. 

71. Cusimano N, Renner SS: Slowdowns in diversification rates from real 
phytogenies may not be real. Sysf Biol 2010, 59:458-464. 

72. Ricklefs RE: Historical and ecological dimensions of global patterns in 
plant diversity. Biol Skr 2005, 55:583-603. 

73. Hoorn C, Wesselingh FP, ter Steege H, Bermudez MA, Mora A, Sevink J, 
Sanmarrjn I, Sanchez-Meseguer A, Anderson CL, Figueiredo JP, Jaramillo C, 
Riff D, Negri FR, Hooghiemstra H, Lundberg J, Stadler T, Sarkinen T, 
Antonelli A: Amazonia through time: Andean uplift, climate change, 
landscape evolution, and biodiversity. Science 2010, 330:927-931. 

74. Couvreur TLP, Chatrou LW, Sosef MSM, Richardson JE: Molecular 
phylogenetics reveal multiple tertiary vicariance origins of the African 
rain forest trees. BMC Biol 2008, 6:54. 

75. Schulte P, Alegret L, Arenillas I, Arz JA, Barton PJ, Bown PR, Bralower TJ, 
Christeson GL, Claeys P, Cockell CS, Collins GS, Deutsch A, Goldin TJ, 
Goto K, Grajales-Nishimura JM, Grieve RA, Gulick SP, Johnson KR, 
Kiessling W, Koeberl C, Kring DA, MacLeod KG, Matsui T, Melosh J, 
Montanari A, Morgan JV, Neal CR, Nichols DJ, Norris RD, Pierazzo E, et al: 
The Chicxulub asteroid impact and mass extinction at the Cretaceous- 
Paleogene boundary. Science 2010, 327:1214-1218. 

76. Wilson R, Heinrichs J, Hentschel J, Gradstein SR, Schneider H: Steady 
diversification of derived liverworts under tertiary climatic fluctuations. 
Biol Lett 2007, 3:566-569. 

77. Venditti C, Meade A, Pagel M: Phytogenies reveal new interpretation of 
speciation and the Red Queen. Nature 2010, 463:349-352. 

78. Morion H, Potts M, Plotkin J: Inferring the dynamics of diversification: a 
coalescent approach. PLoS Biol 2010, 8:el 000493. 

79. McKenna DD, Farrell BD: Tropical forests are both evolutionary cradles 
and museums of leaf beetle diversity. Proc Natl Acad Sci USA 2006, 
103:10947-10951. 

80. Dransfield J, Uhl NW, Asmussen CB, Baker JW, Harley MM, Lewis CE: A new 
phylogenetic classification of the palm family, Arecaceae. Kew Bull 2005, 
60:559-569. 

81. Dransfield J, Rakotoarinivo M, Baker WJ, Bayton RP, Fisher JB, Horn JW, 
Leroy B, Metz X: A new Coryphoid palm genus from Madagascar. Bof J 
Linn Soc 2008, 156:79-91. 

82. Harley MM: A summary of fossil records for Arecaceae. Bof J Linn Soc 
2006, 151:39-67. 

83. Harley MM, Baker WJ: Pollen aperture morphology in Arecaceae: 
application within phylogenetic analyses, and a summary of the fossil 
record of palm-like pollen. Grana 2001, 40:45-77. 

84. Magallon SA: Dating lineages: molecular and paleontological approaches 
to the temporal framework of clades. Int J Plant Sci 2004, 165:S7-S21. 

85. Palmer AR, Geissman J: Geologic Time Scale Boulder, CO: The Geological 
Society of America; 1 999. 

86. Berry EW: In The Upper Cretaceous and Eocene floras of South California and 
Georgia. Volume 84. Reston, VA: United States Geological Survey; 1914. 

87. Pan AD, Jacobs BF, Dransfield J, Baker WJ: The fossil history of palms 
(Arecaceae) in Africa and new records from the late Oligocene (28-27 
Mya) of North-Western Ethiopia. Bot J Linn Soc 2005, 151:69-81. 

88. Schrank E: Palynology of the Yesomma Formation in Northern Somalia: a 
study of pollen spores and associated phytoplankton from the late 
Cretaceous Palmae Province. Palaeontogr Abt B Palaeophytol 1994, 
231:63-112. 

89. Rull V: Biogeographical and evolutionary considerations of Mauritia 
(Arecaceae), based on palynological evidence. Rev Paleobot Palyno 1998, 
100:109-122. 

90. Campbell JD, Fordyce RE, Grebneff A, Maxwell PA: Fossil coconuts from 
mid-Cenozoic shallow marine sediments in southern New Zealand. 
Geo/ Soc New Zeal Misc Publ 2000, 108A:21 . 

91. Medus J: Palynology of the Tertiary sediments from meridional Senegal 
[translated from French]. Pollen et Spores 1975, 17:545-601. 

92. Gomez-Navarro C, Jaramillo C, Herrera F, Wing SL, Callejas R: Palms 
(Arecaceae) from a Paleocene rainforest of northern Colombia. Amer J 
Bot 2009, 96:1300-1312. 



Couvreur ef al. BMC Biology 201 1, 9:44 
http://www.biomedcentral.eom/1 741 -7007/9/44 



Page 12 of 12 



93. Gee CT: On the fossil occurrence of the mangrove palm Nypa. In 

Paleofloristics and Paleoclimatic Changes in the Cretaceous and Tertiary. 
Edited by: Knobloch E, Kvacek Z. Prague, Czech Republic: Geological Survey 
Publisher: 1990:315-319. 

94. Ediger VS, Bati Z, Alisan C: Paleopalynology and paleoecology of Calamus- 
like disulcate pollen grains. Rev Paleobot Paiyno 1990, 62:97-105. 

95. Baker WJ, Dransfield J: Towards a biogeographic explanation of the 
calamoid palms. In Monocots: Systematics and Evolution. Edited by: Wilson 
KL, Morrison DA. Melbourne, Australia: CSIRO; 2000:545-553. 

96. Muller J: Palynology of the Pedawan and Plateau sandstone formations 
(Cretaceous - Eocene) in Sarawak, Malaysia. Micropaleontology 1968, 
14:1-37. 

97. Friis EM, Pedersen KR, Crane PR: Araceae from the Early Cretaceous of 
Portugal: evidence on the emergence of monocotyledons. Proc Natl Acad 
Sci USA 2004, 101:16565-16570. 

98. Drummond AJ, Rambaut A: BEAST: bayesian evolutionary analysis by 
sampling trees. BMC Evol Biol 2007, 7:214. 

99. Rambaut A, Drummond AJ: Tracer, version 1.4.[http://evolve.zoo.ox.ac.uk/]. 

100. Suchard MA, Weiss RE, Sinsheimer JS: Bayesian selection of continuous- 
time Markov chain evolutionary models. Mol Biol Evol 2001, 18:1001-1013. 

101. Ricklefs RE: Phylogenetic perspectives on patterns of regional and local 
species richness. In Tropical Rainforests: Past, Present, & Future. Edited by: 
Bermingham E, Dick CW, Moritz C. Chicago, IL: University of Chicago Press; 
2005:16-40. 

102. Rabosky DL: Likelihood methods for detecting temporal shifts in 
diversification rates. Evolution 2006, 60:1152-1164. 

103. Rabosky DL: LASER: a maximum likelihood toolkit for detecting temporal 
shits in diversification rates from molecular phylogenies. Evol Bioinform 
Online 2006, 2:257-260. 

104. Yule GU: In A mathematical theory of evolution based on the conclusions of 
Dr. }.C. Willis. Volume 213. Philos Trans R Soc Lond, B; 1924:21-87. 

105. Govaerts R, Dransfield J: World checklist of palms Kew, UK: Royal Botanic 
Gardens: 2005. 

106. Tiffney BH, Manchester SR: The use of geological and paleontological 
evidence in evaluating plant phylogeographic hypotheses in the 
Northern Hemisphere tertiary, frit J Plant Sci 2001, 162:S3-S17. 

107. Ali JR, Aitchison JC: Gondwana to Asia: plate tectonics, paleogeography 
and the biological connectivity of the Indian sub-continent from the 
Middle Jurassic through latest Eocene (166-35 Ma). Earth Sci Rev 2008, 
88:145-166. 

108. Morley RJ: Interplate dispersal paths for megathermal angiosperms. 

Perspect Plant Ecol Evol Syst 2003, 6:5-20. 

109. Ronquist F: DIVA (Dispersal Vicariance Analysis) l.i User's Manual 1996 
[http://ib.berkeley.edu/courses/ib200b/labs/labl5/dlvawin/DMANUAL_v2. 
DOC]. 

110. Meerow AW, Noblick L, Borrone JW, Couvreur TLP, Mauro-Herrera M, 
Hahn WJ, Kuhn DN, Nakamura K, Oleas NH, Schnell RJ: Phylogenetic 
analysis of seven WRKY genes across the Palm subtribe Attaleinae 
(Areceaceae) identifies Syagrus as sister group of the Coconut. PLoS ONE 
2009, 4:e7353. 

111. Maddison WP, Maddison DR: Mesquite: A modular system for 
evolutionary analysis, version 2.7. [http://mesquiteproject.org]. 

112. Olson DM, Dinersteln E, Wikramanayake ED, Burgess ND, Powell GVN, 
Underwood EC, D'Amico JA, Itoua I, Strand HE, Morrison JC, Loucks C, 
Allnutt T, Ricketts TH, Kura Y, Wettengel W, Kassem K: Terrestrial 
ecoregions of the world: a new map of life on Earth. Bioscience 2001, 
51:933-938. 



doi: 1 0. 1 1 86/1 741-7007-9-44 

Cite this article as: Couvreur ef ai: Origin and global diversification 
patterns of tropical rain forests: inferences from a complete genus-level 
phylogeny of palms. BMC Biology 201 1 9:44. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at f \ RioM _-| rpntr ,i 

www.biomedcentral.com/submit \ J Blomea central 



